Heavy quark free energies for three quark systems at finite temperature 
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We study the free energy of static three quark systems in singlet, octet, decuplet and average 
color channels in the quenched approximation and in 2-fiavor QCD at finite temperature. We show 
that in the high temperature phase singlet and decuplet free energies of three quark systems are 
well described by the sum of the free energies of three diquark systems plus self energy contributions 
of the three quarks. In the confining low temperature phase we find evidence for a Y-shaped fiux 
tube in SU(3) pure gauge theory, which is less evident in 2-flavor QCD due to the onset of string 
breaking. We also compare the short distance behavior of octet and decuplet free energies to the 
free energies of single static quarks in the corresponding color representations. 

PACS numbers; 



I. INTRODUCTION 

The study of baryonic systems composed of three static 
quarks sheds light on the internal structure of the baryon. 
Lattice calculations have been carried out addressing this 
question for quite some time now [J, H, [3, ll, ^, 6., /7|, . Of 
particular interest is the question whether a genuine three 
body force exists between the quarks below the critical 
temperature and how the system behaves at finite tem- 
perature. At zero temperature the Y-string shape of the 
color flux tube in baryonic systems is supported by recent 
calculations in lattice QCD [1, d, Q . At finite temper- 
ature work so far has concentrated on simulations using 
maximal abelian gauge [1, [Ti] , showing the existence of a 
Y-shaped string as well. The Y-Ansatz is also supported 
by results obtainde with the field correlator method 0. 

Recently, there have been speculations that hadronic 
bound states, e. g. also colored QQQ-sia,ies, may still ex- 
ist in the high temperature deconfined phase of QCD and 
that they could be responsible for the observed strong 
interactions in the QGP near Tc This makes it in- 
teresting to analyze the force in three quark systems also 
in non-singiet color configurations that might exist in an 
overall color neutral thermal medium. 

In this paper, we examine the question of the flux tube 
shape of the QQQ-singlet state below the critical temper- 
ature in the quenched approximation of QCD and discuss 
the appearance of string breaking in 2-flavor QCD. More- 
over, above Tc we explore the relation between the free 
energies in different color channels of the baryonic sys- 
tem and the free energies of QQ-subsystems, and analyze 
the screening of octet and decuplet free energies at small 
distances. 

This paper is organised as follows: In section [ll] we 
present the basic observables we use for calculations in 
different color channels of the heavy three quark system 
in terms of correlation functions of thermal Wilson lines. 
In section Hill we discuss the perturbative behavior of free 
energies of three quark systems at short distances. In sec- 



tion IIVI we present the details of our simulation, includ- 
ing calculation techniques, gauge fixing and renormalisa- 
tion procedures. In section |V] we discuss the behavior of 
QQQ-bee energies obtained in SU(3) pure gauge theory 
on 32'^ X 4, 8 lattices and compare results above Tc with 
free energies of QQ-systems. Furthermore, we examine 
the string shape below Tc- Section IVTl is devoted to a dis- 
cussion of our results in 2-flavor QCD. In section Ivnl we 
analyze the screening of octet and decuplet free energies 
at small distances. Finally, in section IVIIIl we conclude. 



II. STATIC THREE QUARK SYSTEMS IN 
DIFFERENT COLOR REPRESENTATIONS 

In the following we will construct observables for the 
free energies of static three quark systems in different 
color channels and introduce the corresponding relation 
to the expectation values of correlation functions of ther- 
mal Wilson lines. Since we are interested in QCD, we 
restrict ourselves to examine the SU(3) gauge group. 

The state of a three quark system as the product of 
irreducable representations of three quarks in color space 
can be decomposed into symmetry states 

303(^3= 108® 8'e 10. (1) 

The singlet is totally anti-symmetric, the first octet anti- 
symmetric in the first and second, the second octet in the 
second and third component and the decuplet is totally 
symmetric. 

We now construct observables for static three quark 
systems represented by suitable combinations of Wilson 
lines. The derivation is similiar to that for two quark sys- 
tems but more elaborate. We start by defining 
the thermal Wilson line 

Nr-l 

L(x) =. Yl Ui{^,Xi). (2) 
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The Polyakov loop is then obtained by Z(x) = TrL(x), 
where the trace is normahsed such that Tr 1 = 3. For 
the three quark correlation function of thermal Wilson 
lines "'^ Li, L2, -L3, we have 



L1L2L3 — C, 



(xi,X2,X3)P3 



(3) 



We find 



(27 Tr LiTr iaTr L3 - 9 Tr LiTr {L2L3) 

- 9 Tr iaTr {L1L3) - 9 Tr LgTr {L1L2) 

+ 3Tt{LiL2L3)+3Ti-{LiL3L2)) (10) 



where s G {1,8, 8', 10} stands for the symmetry states, 
Ps for the projectors on these states and Cqqq for the 
three point correlator of the thermal Wilson line belong- 
ing to that symmetry state. Here we have suppressed the 
dependence on the temperature T . 

Denoting the components of the thermal Wilson lines 
Lf,L2™,L^"', where i,j,k,l,m,n — 1,2,3, we find the 
projectors by using the Fierz-identity 



QQQ 



= — (27 Tr LiTr LaTr L3 + 9 Tr LiTr {L2L3) 

-9TrL3Tr(Lii2) -3Tr(LiL3^2)) (11) 



QQQ 



= — (27 Tr LiTr L2Tr L3 + 9 Tr L3Tr {L1L2) 

-9TrLiTr(L2i3) -3Tr(LiL2L3)) (12) 



(4) 
(5) 
(6) 
(7) 



where the sums are over all permutations a of the indices. 
The Ps satisfy the usual relations for projectors 



^10 
^QQQ 



— (27 Tr LiTr L2Tr L3 + 9 Tr LiTr (L2i3) 

+ 9 Tr L2Tr {L1L3) + 9 Tr PgTr {L1L2) 
+ 3 Tr {L1L2L3) + 3 Tr {L1L3L2)) . 



(13) 



Finally we obtain for the color averaged correlator of the 
three quark system, Cnnn. the relation 



'-'QQQ 



= — Tr L^ Tr PoTr L3. 
27 



+ 27^QQQ 



10^10 
2'j'^QQQ 

(14) 



We note that the two octet correlators Cqqq and Cqqq 
are the same when calculated on the lattice. The free 
energy of the symmetry state s can be obtained from the 
correlator in the usual way 



P. 



E 



Ps^l und P.Pt = 



(8) 



ioT s ^ t and s,t € {1,8, 8', 10}. The desired three quark 
correlators in the different color channels are now ob- 
tained by applying 



FQQQiT) = -Tlr,{C^QQiT)), 



(15) 



where (•) stands for thermal averages taken after gauge 
fixing. We have suppressed the position dependence on 
both sides of (fT5)) . 

In the following section we will discuss briefly the be- 
havior of the free energies of the three quark system for 
small couplings and short distances. 



c 



QQQ 



Tr ( PsLiL2L, 
Ti-Ps 



(9) 



^ For convenience we write Li instead of L(xi). 



III. 



f'qqq at short distances 



In the perturbative expansion of the free energy of 
QQQ-systems the contribution of the three gluon vertex 
vanishes for symmetry reasons [isj . Therefore, neglect- 
ing self energy contributions, to order the free energy 
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Fqqq decomposes into the sum of three diquark free en- 
ergies Fqq, which can be in an anti-symmetric antitriplet 

{t — 3) or in a symmetric sextet {t = 6) state. The be- 
havior of these diquark free energies at small distances 
to lowest order are Coulombic 



F\R,T) = C,{t)-, 



(16) 



where a = |^ and R is the separation of the static quarks 
in the QQ-system. The Casimir factor C2 = Tr (^1*2) de- 
pends on the color channel t the diquark system is in and 
can be found in tabic [H In order to obey the permuta- 
tion relations given for the QQQ-system in ([4]) - ([7]) , the 
QQQ-singlet state must be composed of QQ-anti-triplets, 
the QQQ-decuplet of QQ-sextets and the QQQ-octets 
are a mixture of QQ-anti-triplets and -sextets. Thus, we 
have for small distances 



i<j 



fci(T) 



QQQ 



F 



QQQ 



(R,T) 



(R,T) 



2 a 
1 a ^ 

3 i?i2 



1 a 
1 a 



_^ 1 a 
3 i?23 

2 a 

3 i?23 



F, 
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QQQ 



(R,T) = Y^l^ + koiT), 



(17) 

(18) 

+ ~k8'{T) 

(19) 
(20) 



where Rij denotes the distance between the ith and jth 
quark (see sec. lIVBl for details), R = (i?i2, ^13, -^23)- 
The self energy contributions ki with i = 1,8, 8', 10 
are temperature dependent and at small separations re- 
lated to the free energy of a single static quark in the 
corresponding color state, Fq \t) (see sec. IVIII for de- 
tails). Thus fci(r) = 0, ks{T) = h'iT) = F^\t) 

and fcio(T) = Fq^\t). As can be seen easily from 
(IT7|) and PU)1 . the singlet state of the QQQ-system is 
attractive, because the QQ-antitriplet is, and the decu- 
plet state is repulsive, as the QQ-sextet is. Moreover, 
the QQQ-singlet free energy is temperature independent 
at small distances like the QQ-singlet free energy is [l5| . 
For the octet channels the situation is more complicated 
due to the contribution of both anti-triplet and sextet 
free energies. We can, however, compute the lowest or- 
der behavior of the QQQ free energies for equilateral ge- 
ometries for all color channels. Table U summarizes the 
Casimir factors for the free energies of the different color 
channels for all the quark systems QQ, QQ and QQQ. 
One recovers the lowest order behavior of the free energy 
in some symmetry state t (equilateral geomtries for the 
QQQ-systems only) that depends on the separation of 
the static quarks by using (|16p . where R is now the sep- 
aration of the static quarks in the QQ- and QQ-systems 



system 1 



6 



8 10 



QQ -4/3 +1/6 

QQ -2/3 +1/3 

QQQ^ -2 -1/2* +1 



TABLE I: Casimir factor for color symmetry state t, C2{t) 
Tr {T?-Ti). 

f equilateral geometries, * average of both octets. 



and the edge length of the equilateral geometry in the 
QQQ case. See sec. IIVBI for details on the geometric 
configuration of the three static quarks. We see, that the 
average of both QQQ-octet free energies (which is acces- 
sible to lattice calculations) for equilateral geometries is 
expected to be weakly attractive. 

At larger distances, the pairwise interactions of the 
static quarks should receive temperature dependent con- 
tributions and the Coulomb terms in ([TT]) - ([201) have to 
be replaced by the full diquark free energies Fqq (i?, T) 
of the corresponding color channel t. In this case, pT)) - 
(PO)) generalize to 



F^QQ(R,r) = ^F3Q(i?,,-,T) + fci(r) 



i<3 



+J^^q(^23) + A:8(T) 
Pqqq(^^T) = -^^00(^12) - 7-^Qq(^13) 



{R23) + ks'iT) 



+Fqq 



i<j 



F^Q{R,,,T) + k,o{T), 



(21) 



(22) 



(23) 



(24) 



where the ki{T) with i — 1,8, 8', 10 account for the 
self energy contributions not included in the diquark 
free energies. It was shown in [l3| that the residual T- 
dependence of the QQ-anti-triplet free energy at small 
distances in the deconfined phase can be removed by 
substracting the free energy of a single static quark 
Fq{T) = ^limR^oo FI^q{R,T). We therefore expect 
fci(r) = — 3Fq(T), which guarantees the cancelation 
of self energy contributions stemming from Fqq, leav- 
ing Fqqq tempcrature independent at small separations. 
Moreover, at large separations, the self energy contribu- 
tions for the QQQ-system should be independent from 
the particular color channel the system is in, leading to 
fci(r) = ks{T) = ks'{T) = fcio(T) - -3Fq{T). We will 
show in sections IVIVII and IVIII that this reasoning is in- 
deed correct for temperatures T > Tc. Indeed, we will 
show that the decomposition of interactions in 3-quark 
systems in terms of 2-quark interactions as suggested by 



FIG. 1: Inter quark distances. F is the Fermat point of the 
triangle. 



FIG. 2: Calculation of the three point correlation function of 
the Polyakov loop. 



(I2I])-(I211) holds for T > Tc. 

We finally note, that the QQQ-octet free energy ob- 
tained from the lattice is the average of Fqqq and Fqqq , 
as will be shown in sec. IIVI Consequently, we expect the 
lattice QQQ-octet free energy for equilateral geometries 
to obey 



-3Fq(T), 



(25) 



where R is the edge length of the equilateral triangle. 



IV. SIMULATION DETAILS 



B. Calculation technique 

We first fix the notation for the geometries of the three 
quark system. In fig. [T]we show three quarks Qi forming 
a triangle and their distances Rij, where i,j — 1,2,3. 
The perimeter of the triangle is then simply given by 



P 



i<j 



(26) 



This is also the total length of a A-shaped string, i. e. a 
string connecting the three quarks along the edges of the 
triangle. Another possible string shape is a Y-shaped 
string, where the flux tube emanetes from each of the 
three quarks and has a junction at the Fermat point F of 
the triangle. The total length of such a Y-shaped string 
is 



We will now present the details of our simulation. 



A. Simulation 

We used gauge field configurations generated on 32-^ x 4 
and 32'^ x 8 lattices in pure gauge theory with the tree 
level-Symanzik improved gauge action at several cou- 
plings above and below the critical coupling [l^ [13] ■ In 
2-flavor staggered QCD we reexamined configurations on 
a 16'^ X 4 lattice for several different couplings at bare 
quark mass m/T = 0.4, where for fermions the p4-action 
and for the gauge fields again the tree level-Symanzik im- 
proved gauge action were used [H, [l^ . The scale was 
set using the string tension, cr, at T — following f^Tl 
and use it to express all dimensionful observablcs. On all 
configurations we calculated the three point correlation 
functions of Polyakov loops (fTU]) - (fT5|) , (|14p in a manner 
explained below. 

The operators in (fTO| - (|13p are not manifestly gauge 
independent and therefore a gauge fixing procedure must 
be applied to our gauge configurations. We use Coulomb 
gauge for our calculations, in which the singlet free ener- 
gies of the QQ-system are related to a gauge independent 
definition in terms of dressed Polyakov lines Note 
that an operator dependence might still persist within 
this definition [23^. 



L = 



i<3 



(27) 



where Aa is the area of the triangle p . The inner angles 
at the vertices of the triangle are assumed to be smaller 
than which is the case for all triangles considered in 
this work. The angle between any two arms of the Y- 
shaped fiux tube is always For equilateral triangles 
we have the simple relation P = \/3L. In this work 
we examine only isoscele triangles, where we set i?i2 = 
i?i3 — Rs to be the equally long edges and R23 = Rb to 
be the length of the base edge. For equilateral triangles 
we have Rij — R for all i,j = 1, 2, 3. 

The three point correlation functions of the Polyakov 
loop are now obtained as follows (see fig. ^ . At the posi- 
tions of the Qi, we calculate the correlation functions (fTO|) 
- (IT5|) . (|14p and compute the average of the correlation 
function for those Qi with the same {Rij}. We obtain a 
new set of {Rij} by holding Qi fixed, whereas the two 
other vertices of the triangle Q2 and Q3 are moved simul- 
tanously one point of the lattice in one direction (here: 
to the left). The base edge Rb connecting these to points 
preserves thereby its length, which is Rb/a = 7i\/2, where 
a denotes the lattice spacing and 1 < n < ^ is an in- 
teger and describes the number of elementary cells the 
procedure starts with. The two other edges have equal 
lengths Rs/a = \Jm?- n^, where m is another integer 
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F(R,T)/o^ 




0.1 0.2 0.3 0.4 0.5 0.6 



FIG. 3: Free energies in different color channels obtained on 
a 32^ X 8 lattice at T /Tc = 6 from simulations in quenched 
QCD. 

which runs between n < m < ^f- for every n. Therefore 
for every n we obtain one equilateral (n = m) and several 
isosceles (n < m < triangles. We start with n = 1 
and repeat the procedure until n — — We apply this 
method in both directions of all three spatial dimensions 
before sweeping Qi over the entire spatial lattice. We will 
denote the configuration averages of the correlation func- 
tions (Uni) - (HSl) and mi) by C'qqq with s = 1,8, 10, av 
from now on. 



C. Renormalization 

We construct all correlation functions using renor- 
malised thermal Wilson lines L''®"(T), which are obtained 
from the bare ones calculated on lattices with temporal 
extend Nt through multiplicative renormalisation, 

U^%T) = {Z{g^))''^ L{g\Nr), (28) 

where Z{g'^) is the multiplicative renormalisation con- 
stant determined in [l5| , g is the bare coupling and N^. is 
the temporal extent of the lattice L is calculated on. The 
renormalization constants only depend on the bare cou- 
pling 24, 25] (and in addition on the bare quark masses 
in full QCD) and furthermore are equivalent at zero and 
finite temperature [1^. The renormalised three point 
correlation function of the thermal Wilson lines, C"''^"(T), 
is then obtained from the bare one, CqQqig^, Nr), by 

CQQQiT) - {Zig')f''^ CQQQig^Nr). (29) 

V. RESULTS IN PURE GAUGE 

In this section we will analyze the behavior of three 
quark free energies in different color channels above T^. 
We compare results with the perturbative expressions ob- 
tained in section IIIII as well as with results obtained for 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



FIG. 4: F^QQiR,T) and 3 {F^QiR,T) ~ Fq{T)) above 
versus R, the edge length R of the equilateral triangles and 
QQ-distance, respectively. 

QQ-systems in different color channels. Data for the QQ 
free energies have been taken from [l3|. Below Tc we 
examine the string shape of the baryonic system and its 
structure in the different color channels. 



A. Color Channels 

In fig. [3] we show the free energies of three quark sys- 
tems in different color channels and the average free en- 
ergy for the QQQ-system for equilateral triangles of edge 
length R calculated on a 32^ x 8 lattice at T/T^ = 6. 
One can see clearly, that the singlet is strongly, the octet 
weaker attractive and the decuplet repulsive in agree- 
ment with the perturbative findings presented in sec. IIIII 
For large R at a given temperature, all free energies in 
the different color channels approach a common value, 
i. e. the three quarks are screened independently of their 
color orientation. The singlet free energy becomes tem- 
perature independent at small distances and coincides 
with the baryonic T = potential, Vqqq{R), (see also 
fig. H]), which is related to the quark-antiquark poten- 
tial at vanishing temperature by the ratio of the differ- 
ent Casimir operators, i. e. Vqqq{R) = ^Vqq{R) for 
R^QCD ^ 1- We obtain similiar results for all other 
temperatures above Tc. We will discuss the screening 
of octet and decuplet free energies at small distances in 
sec. Ivnl 



B. Free energies of equilateral geometries above 

We now compare the free energies of the QQQ-system 
with the free energy of the QQ-system above Tc- In fig. [J] 
we show F^qq{R,T) and 3F^q{R,T) - 3Fq{T) versus 
the edge length R of the equilateral triangles and the QQ 
distance, respectively. According to (PT|) the second term 
is expected to be equal to the QQQ-singlet free energy at 
least at small distances, where genuine three body forces 
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Fqqq(R,T)/o"^ □ 3(F^q(R,T)-Fq(T))/o^'2 





FIG. 5: FlfQQiR,T) and 3 {F^q{R,T) - Fq{T)) above 
versus R, the edge length R of the equilateral triangles and 
QQ-distance, respectively. 



FIG. 7: AF(P,T) = F^%q{P,T) - J:^^^F^q{R,,,T) + 
SFq{T) above Tc versus the perimeter P for all geometries 
calculated. 





FIG. 6: AF{P,T) ^ F^QQ{P,T)^j:^Fl,g{R,j,T)+3FQ{T) 
above Tc versus the perimeter P for all geometries calculated. 



are negligible. For all temperatures above we indeed 
see that F^qq{R,T) and 3F^QiR,T) - 3Fq{T) coin- 
cide throughout the entire distance interval. In fig. [5] we 
show analogously Fqqq(-R, T) and SFj^giR, T)-3Fq(T). 
Again we observe that both observables do coincide. We 
obtain similar plots for Fqqq{R,T) in accordance with 
13. 



C. Free energies of isosceles geometries above Tc 

In order to test whether the simple relation between 
free energies of a 3 quark system and the free energies of a 
2 quark system plus self energy terms also holds for other 
geometries and temperatures above Tc, we calculate the 
difference AF(P, T) = F^qq{P, r)-E,<, F'^q{R,„T)+ 
3Fq. This is shown in fig.Hl If (|2ip also holds for these 
geometries, then AF{P, T) should vanish, which is full- 
filled to a very good degree for all perimeters except the 
smallest ones. Again the same holds true also for the 
QQQ-decuplet free energy as is evident from fig. [T] This 
is possibly due to the effect, that the self energy of the 



FIG. 8: AF{P,T) = FhQQiP.T) - T,.^^ Fi^giR.j ,T) be- 
low Tc versus the perimeter P for all geometries calculated. 
Ai^(Pniin,r) has been set to zero, where Pmin is the smallest 
perimeter calculated. 



system is no longer the sum of the self energy of the in- 
dividual static quarks at short distances but rather that 
of a single static quark in the corresponding color repre- 
sentation (see sec. IVII[) . 

Having now established that (PT|) and hold for all 
isoscele geometries above Tc calculated in this work, it is 
clear that the screening length of the three quark system 
in the singlet state is the same as that of the antitriplet 
diquark state, which itself has been shown to be equal to 
that of the QQ singlet state jlJi] , therefore refiecting the 
properties of the thermal medium rather than that of a 
particular hadronic system. 

Furthermore, we can apply the relations used in to 
calculate entropy and internal energy from free energies 
of a QQ system also for the QQQ system. The entropy 
in the color singlet channel is defined by 



{P,T) 



dT 



(30) 
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FIG. 9: F^qq{.,T) at T/Tc = 0.925 versus the Y-string 
length L and the perimeter P for equilateral geometries (eq) 
and geometries with base length Rh^f^ = n^/a\/2, where n is 
given in the legend. 



and the internal energy in this channel is 



dT 



Wc then obtain for the entropy 



(31) 



(32) 



i. e. the entropy contribution in the free energy of the 
QQQ singlet color channel vanishes at small distances. 
Moreover, we find for the internal energy, that 



(33) 



where Vqqq{P) = 2E»<j ^qq(%) and Yggir) is the 
quark- anti-quark potential at T = 0. This means that 
the internal energy becomes T-independent for small 
perimeters and thus, together with f30|) . the singlet free 
energy itself is T-independent for small perimeters, as 
was already seen in sec. IV Al This behavior is already 
known for QQ singlet free energies and reflects the fact 
that the baryonic system is less and less effected by the 
surrounding thermal medium when going to smaller and 
smaller perimeters. 



D. Free energies below Tc 

We now examine the QQQ free energies below Tc- We 
start by looking at the relation between the QQQ-singlet 
and the QQ-antitriplet free energies. If Fqqq can be 
expressed in terms of the sum of FQQ{Rij ,T) also be- 
low Tc, i. e. if a A-ansatz for the flux tube shape to- 
gether with the same string tension holds, then Fqqq 
is a function of the perimeter only and AF(P, T) — 
Fqqq{P,T) - j:i<,FlQ{R.j,T) should be equal to a 
T-dependent constant k{T) for all P. In fig. [5] we show 



FIG. 10: Free energies of the color channels versus the Y- 
string length L at T/Tc = 0.925 for all geometries calculated. 



AF{P,T)/T for T <Tc versus the perimeter, P, for dif- 
ferent geometries. Here we have set AF{Pnuu,T) — 0, 
where Pmin is the smallest perimeter calculated. We 
can see clear deviations of AF{P, T) to smaller values, 
most strongly for the lowest temperatures, deviations 
from zero becoming smaller with growing temperature. 
For T/Tc = 0.986 wc have AF{P,T) « 0. Hence, 
Fqqq{P,T) can not be expressed in terms of the sum 

of Fqq [Rij , T) except close to the critical temperature. 
We are left with two possibilities now. First, we could 
still have a A-shaped flux tube, but with a different string 
tension than that observed for the QQ-antitriplet. Since 
we see AF{P, T) < 0, we would expect the string tension 
to be smaller than in the QQ-antitriplet. In this case the 
perimeter P would still be the right distance measure 
for the QQQ-singlet free energy, i. e. Fqqq should be 
a smooth function of P for all geometries. Or, secondly, 
the flux tube is Y-shaped and L is the right distance mea- 
sure. In this case, Fqqq should be a smooth function of 
L for all geometries in which a flux tube can form. 

To elaborate more on the shape of the flux tube, we 
take a closer look at the QQQ-singlet free energy at a 
particular temperature below Tc- We analyze the free 
energy for different geometries as a function of L and 
P. For equilateral triangles a simple geometrical relation 
exists between the length of a Y-shaped flux tube L and 
the length of a A-shaped flux tube P, which is P = ^/SL. 
Hence for equilateral geometries the QQQ-singlet free en- 
ergy is a smooth function in both Ansatze. For more gen- 
eral geometries like the isoscele triangles we calculated, 
no simple relation between L and P exists. This may 
help to clarify the situation more directly. Therefore we 
plot the QQQ-singlet free energy at T/Tc — 0.925 versus 
the Y-string length L and the perimeter P. This is shown 
in flg.[9]for equilateral geometries and also for isoscele ge- 
ometries up to n = 5 (see sec. lIVB|) . We observe that the 
QQQ-singlet free energy is indeed a smooth function of 
L for all geometries calculated. On the other hand, when 
the free energy is plotted versus the perimeter P differ- 
ent branches become visible depending on the geometry. 
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FIG. 11: Fqqq{L,T) below Tc versus the Y-string length L FIG. 12: Free energies in different color channels in 2-flavor 
for all geometries. QCD from calculations on a 16^ x 4 lattice at T/Tc = 1.99. 



This is most prominent for the n = 1 triangles. Therefore 
we find strong hints, that the shape of the flux tube in 
the QQQ-singlet system is indeed that of the Y-ansatz 
and its string length L is the right distance measure for 
the system. 

Having estabhshed this, we take a look at the other 
color channels below Tc and analyze whether they display 
a smooth behavior as function of L as weh. In fig. [TUlwe 
show the different cofor channels of the QQQ-bee energy 
at T/Tc — 0.925 over the Y-string length L. We see again 
that the singlet is the most attractive channel followed by 
the average free energy, which is also a smooth function 
of L. The octet channel is still attractive, but weaker so 
than the average free energy. The decuplet free energy is 
attractive for large L but becomes flat at smaller L, pos- 
sibly hinting at a turnover and at a repulsive behavior at 
even smaller L. Both the decuplet and the octet channel 
are not smooth functions over either distances L and P, 
but become volatile for Ly/a>2. This suggests that no 
flux tube forms in these two channels, besides that the 
octet free energy shows an overall attractive behavior. 

Finally, we examine more closely the temperature de- 
pendence of Fqqq . In flg. [11] we show the renormalised 
QQQ-singlet free energy versus L for all temperatures 
below Tc. We observe F^qq{L,T) to nearly coincide 
at distances Ly/a<4: for all temperatures T/Tc < 0.959. 
We see that at large distances the free energies approach 
F{L, T) = aL represented by the black line in flg. [Til 
where a denotes the string tension in a QQ-system. The 
phenomenon of an almost T-independcnt string tension 
has also been observed in Q, using a flt ansatz in full 
QCD. As was aheady seen in flg. [3 at T/Tc = 0.986 we 
observe deviations to smaller values. This might be due 
to flux tube broadening setting in near Tc, leading to an 
overlapping of the three branches of the Y-shaped flux 
tube. Filling the space between the three static quarks 
the separation dependence of the singlet free energy is 
effectively described by the A- Ansatz at these tempera- 
tures. 



VI. RESULTS IN 2-FLAVOR QCD 

In this section we analyse the QQQ-free energies ob- 
tained on a 16^ X 4 lattice in 2-flavor QCD with a bare 
quark mass of ^ = 0.4 . For further details, see se c. IIVI 
Data for the QQ free energies is again taken from [lj|. 

A. Color Channels 

In flg. [13 we show the free energies of three quark sys- 
tems with equilateral geometries in different color chan- 
nels at T/Tc = 1.99. Like in the pure gauge case, we 
observe the singlet to be strongly, the octet weaker at- 
tractive and the decuplet to be repulsive. For large dis- 
tances R the free energies in all color channels approach 
a common value. Again, we obtain similar results for all 
T>Tc. 

B. Singlet free energy 

We now compare the free energy of the QQQ-singlet 
with the free energy of the QQ-antitriplet for 2-flavor 
QCD. In flg. [131 we show F^qq{R,T) and 3F^q{R,T) - 
3Fq(T) versus the edge length R of the equilateral tri- 
angles and QQ distance, respectively. As in pure gauge 
theory, both curves coincide for all temperatures above 
Tc. Thus we flnd the QQQ-singlet free energy of equi- 
lateral geometries can be described as the sum of three 
QQ-antitriplet free energies plus self energy terms as well 
in 2-flavor QCD above Tc. 

Below Tc we expect to see string breaking also in the 
QQQ-singlet free energy, the breaking mechanism being 
more involved than in the QQ case [T]. Indeed, except 
at T/Tc — 0.97, where both quantities agree, we see de- 
viations for the QQQ-singlet from the simple relation to 
the QQ free energies. In fig. [14] we compare Fqqq{L,T) 
in 2-flavor QCD (T/Tc = 0.88, 0.91) and pure gauge the- 
ory (T/Tc = 0.874, 0.907) for aU geometries calculated. 
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FIG. 13: F^qq{R,T) and 3{F'^q{R,T) - Fq{T)) versus R, 
the edge length ol the equilateral triangles and QQ-distance, 
respectively. 



FIG. 15: F^%q{R,T) and 3{F^q{R,T) - Fq{T)) versus /?, 
the edge length of the equilateral triangles and QQ-distance, 
respectively. 




1 1.5 2 2.5 3 3.5 4 4.5 5 



FIG. 14: Comparison of F^qq{L,T) in 2-flavor QCD (open 
symbols) and SU(3) pure gauge theory (filled symbols) for 
two temperatures below Tc. 



geometries in fig. 1151 Like in pure gauge theory, we 
observe both quantities to coincide for all temperatures 
above Tc and distances R. The situation below Tc is 
analogous to the singlet case. At the temperature closest 
to Tc both quantities coincide for all R, at smaller T we 
observe deviations to smaller values for Fqqq {R, T) . 

Hence we see, that the singlet (decuplet) QQQ free 
energy can be described as a sum of antitriplet (sextet) 
QQ free energies and the self energy contributions of the 
three quarks above Tc (Eqs (pTjl and (|24p ) also in 2-flavor 
QCD. Again, we obtain similar plots for Fqqq{R,T) in 
accordance with (|^. 



Here the pure gauge data have been obtained from sim- 
ulations on a 32^ x 4 lattice. The Nf — 2 free energies 
start to deviate from the pure gauge result already at 
short distances to smaller values and eventually become 
flat. Specifying a definite value for the string breaking 
distance is quite difficult given the present data. Never- 
theless we can give a rough estimate for the distance at 
which the pure gauge free energies assumes the asymp- 
totic value of the 2-flavor free energies, which is Ly/a sa 3. 
For higher T, this distance becomes smaller. The string 
breaking distance for the QQ singlet free energies has 
been determined in as the distance where the T = 
potential assumes the large distance asymptotic value of 
the free energy. This lead to values between 2.8/y/a for 
T/Tc = 0.874 and 1.7/Vct close to Tc. 



C. Decuplet free energies 

For the QQQ-decuplet free energies we show 
F^qq{R,T) and 3{F^q{R,T) - Fq{T)) for equilateral 



VII. SCREENING OF OCTET AND DECUPLET 
FREE ENERGIES 

As was shown in [l^, the screening of diquark free 
energies in the anti-triplet and sextet representation at 
small distances becomes identical to the screening of a 
single fermion in a color anti-triplet and sextet represen- 
tation, respectively. This holds for all temperatures in 
2-flavor QCD and in the deconflnement phase of SU(3) 
pure gauge theory. 

We will now show, that analogous statements are true 
in the deconflnement phase for octet free energies of the 
QQ-system and for octet and decuplet free energies of 
the QQQ-system, as we have proposed in sec. IHII 

When combining singlet and octet free energies or sin- 
glet and decuplet free energies in proportion to the corre- 
sponding Casimir factors (cf. Tab.|T| such that the repul- 
sive and attractive Coulombic contributions cancel, we 
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FIG. 16: Comparison of singlet and octet free energies at 
small distances with the free energy of a static fermion in 
the adjoint representation. We have used F^-^{R = a,T) = 

F^^iR = a,T) +8F^^{R = a,T) and i^QQQ(i?mi„, T) = 
-2F^QQ{Rn.in,T)+8F^QQiRmin,T). This particular combi- 
nations of QQ- and QQQ-iree energies shown eliminate the 
short distance Coulomb terms contributing to them according 
to p4[) - p6[) . Open symbols: Nt = 4 data, closed symbols: 
A''t- — 8 data. For the QQQ-iree energies we used equilateral 
geometries with edge length Tlmin = 



expect to find 

lim (i^(l^(^,r) + 8F(g(i?,T) 
i?.^o V QQ QQ 



(34) 
(35) 
(36) 



and 



1™ Kqq(^'^) + ^PQQQiR^T)) = 2^;(^")(r) 

(37) 

where R stands for the edge length of the equilateral ge- 
ometries in the case of QQQ free energies; Fq\t) and 

Fq^^ (T) are the free energies of a static quark in a color 
octet and decuplet representation, respectively. The lat- 
ter are given by the corresponding Polyakov loop expec- 
tation values, 



where D = 8, 10 and 

L8(x) = i(9|i(x)p-l) 



(38) 



(39) 



Lio(x) = -(27L(x)3-18|L(x)|2 + l) (40) 



FIG. 17: Comparison of singlet and decuplet free energies at 
small distances with the free energy of a static fermion in 
the decuplet representation. We have used ^qqq (-Rmin, T) = 
FQQQ{Rmin,T) + 2F^°QQ(R^in,T) . This particular combina- 
tion of QQQ-iree energies shown eliminates the short distance 
Coulomb terms contributing to them accroding to (|37p . For 
the QQQ-iree energies we used Nt = 4 data and again equi- 
lateral geometries with edge length -Rmin = V2a. 



are obtained through group theoretical relations. The 
Polyakov loops (ig) ^^'^ (^lo) have been computed in 
[131 ■ We note here, that these quantities - unlike the 
Polyakov loop in the fundamental representation - do 
not have to vanish in the confinement phase of SU(3) 
pure gauge theory due to their vanishing triality. 

(HH) - ((57|) are indeed fulfilled for all temperatures 
above Tc as can be seen from Figs. fTHl and fTTl In partic- 
ular, the combinations of singlet and octet free energies 
of the QQ- and the QQQ-system coincide in the decon- 
finement phase (|55|) , thus showing the screening at small 
distances to depend only on the color representation and 
not on the particular static quark system under consider- 
ation. Below Tc deviations become apparent in Fig. [161 
The reason probably is, that the minimal distances used 
in the present analysis {R — a for the QQ-system and 
Rmin — V2a for the QQQ-system respectively) are still 
too large at small temperatures to be a good approxima- 
tion for the small distance limit in ([M)) - ([57)) . In fact, the 
Nr — 8 data, where the lattice spacing a is only half as 
large at the same value of the temperature, agree much 
better at least close to Tr. 



VIII. CONCLUSIONS 

We have calculated the free energy of static three quark 
systems in different color channels in the quenched ap- 
proximation and in 2-flavor QCD at finite temperature. 
We have shown that above the critical temperature the 
singlet and decuplet free energies of the three quark sys- 
tem can be described by the sum of the free energies of 
the corresponding diquark system plus self energy con- 
tributions. Therefore the screening of the singlet QQQ 
free energies is comparable to that of the QQ free en- 
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orgies reflecting the screening properties of the thermal 
medium. Below T^. we found evidence for a Y-shaped flux 
tube in the quenched approximation for the singlet and 
color average channel. The string tension agrees with 
that deduced from QQ-systems at T = and has been 
found to be independent of temperature for T/Tc<0.959 
in the singlet color channel. In 2-flavor QCD we observed 
string breaking in the QQQ singlet free energies at dis- 
tances around 3/\/o"- Moreover, we were able to show 
that at short distances the screening properties of octet 
free energies of the QQ- and the QQQ-system and the 
decuplet free energies of the QQQ-system coincide with 



those of a single static quark in the corresponding color 
state above T^. 
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